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Abstract 

This  paper  presents  a  neural  networks  (NN)-based  active  and  reactive  power  controller  of  a  stand-alone  proton  exchange  membrane 
(PEM)  fuel  cell  power  plant  (FCPP).  The  controller  modifies  the  inverter  modulation  index  and  the  phase  angle  of  the  ac  output  voltage  to 
control  the  active  and  reactive  power  output  from  an  FCPP  to  match  the  terminal  load.  The  control  actions  are  based  on  feedback  signals 
from  the  terminal  load,  output  voltage  and  hydrogen  input.  The  validity  of  the  controller  is  verified  when  the  FCPP  model  is  used  in 
conjunction  with  the  NN  controller  to  predict  the  response  of  the  power  plant  to:  (a)  computer- simulated  step  changes  in  the  load  active 
and  reactive  power  demand,  and  (b)  actual  active  and  reactive  load  demand  of  a  single  family  residence.  The  response  curves  indicate  the 
load  following  characteristics  of  the  model  and  the  predicted  changes  in  the  analytical  parameters  highlighted  by  the  analysis. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  fuel  cell,  as  a  renewable  energy  source,  is  considered 
one  of  the  most  promising  sources  of  electric  power.  Fuel 
cells  are  not  only  characterized  by  higher  efficiency  than 
conventional  power  plants,  but  they  are  also  environmentally 
clean,  have  extremely  low  emission  of  oxides  of  nitrogen 
and  sulfur,  and  have  very  low  noise. 

The  main  components  of  a  proton  exchange  membrane 
(PEM)  fuel  cell  system  include  the  fuel  processing  unit  or  the 
reformer,  the  fuel  cell  stack  and  the  power  conditioning  unit. 
In  general,  hydrogen  fuel  is  produced  by  processing  some 
hydrocarbon  fuel  such  as  propane,  natural  gas,  or  methanol 
in  the  reformer.  During  fuel  processing,  carbon  monoxide 
is  produced.  Reduction  of  carbon  monoxide  to  acceptable 
levels  is  achieved  by  the  water-gas  shift  reaction. 

The  output  from  the  fuel  cell  is  dc  power.  When  a  fuel  cell 
power  plant  (FCPP)  provides  power  to  a  residential  load,  or 
to  the  electrical  grid,  a  power  conditioning  unit  is  needed. 

The  power  conditioning  unit  is  simply  a  dc/dc  converter 
used  to  raise  the  dc  output  voltage,  which  is  generally  the 
dc  bus  voltage,  followed  by  a  single-phase  or  three-phase 
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dc/ac  inverter.  In  grid  parallel  operation  mode  a  transformer 
is  needed  because  of  the  voltage  difference  between  the  fuel 
cell  power  plant  output  and  the  grid.  Due  to  low  working 
temperature  (80-100  °C)  and  fast  start  up,  PEM  fuel  cells 
are  best  suited  for  residential  and  vehicular  applications. 
The  reasons  behind  selecting  the  PEM  FCPPs  in  this  pa¬ 
per  are:  (a)  fuel  cell  user’s  group  publications  indicate  that 
for  small  residential  application,  where  major  heating  appli¬ 
ances  are  natural  gas  operated,  a  5  kW  unit  is  adequate;  (b) 
a  5  kW  commercial  PEM  FCPP  is  currently  operational  in 
the  authors’  lab. 

Many  models  have  been  proposed  to  simulate  fuel  cells 
in  the  literature.  The  basis  of  a  model  can  be  fluid  dynamics, 
electrochemical  reaction,  heat  transfer  and  thermal  [1-10]. 
To  study  the  transient  response  of  a  grid-independent  PEM 
fuel  cell  power  plant,  this  paper  employs  an  electrochemical 
model  for  a  5  kW  fuel  cell,  which  incorporates  an  external 
reformer  to  generate  hydrogen  from  methane. 

To  enable  the  FCPP  to  conform  to  the  load  changes,  it 
is  essential  to  control  the  active  and  reactive  power  output 
of  the  power  plant.  An  active  and  reactive  power  control 
scheme  has  been  developed  in  Ref.  [11].  The  control  scheme 
is  based  on  controlling  the  ac  output  voltage  phase  angle  to 
control  the  power  flow  from  FCPP  to  the  connected  load.  The 
reactive  power  is  controlled  through  controlling  the  inverter 
modulation  index.  This  paper  uses  the  same  controlling 
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scheme  to  develop  a  neural  network  (NN)-based  controller 
to  control  the  active  and  reactive  power  output  of  the  FCPP. 

The  paper  is  organized  as  follows:  Section  2  introduces 
a  model  for  a  fuel  cell  system.  Section  3  presents  the  NN 
design  and  training  process.  Load  following  test  results  of 
this  model  are  presented  in  Section  4.  Section  5  presents  the 
conclusions. 


Fig.  2.  Reformer  and  the  reformer  controller  model. 


2.  Fuel  cell  system  model 

2.7.  Fuel  cell  model 

In  Ref.  [6],  Padulles  et  al.  introduced  a  model  for  the 
SOFC.  The  model  has  been  modified  to  simulate  a  PEM  fuel 
cell  [7].  This  model  is  based  on  simulating  the  relationship 
between  output  voltage  and  partial  pressures  of  hydrogen, 
oxygen,  and  water.  A  detailed  model  of  the  PEM  fuel  cell 
is  shown  in  Fig.  1.  The  model  parameters  are  as  follows: 

B ,  C  constants  to  simulate  the  activation  overvoltage 

in  PEM  FC  (A-1)  and  (V)  [7,10] 

E0  Standard  no  load  voltage  (V) 

F  Faraday’s  constant  (C/kmol) 

I  stack  current  (A) 

Kr  modeling  constant  (kmol/s  A) 

Kft2  hydrogen  valve  molar  constant  (kmol/atm  s) 

Kk2o  water  valve  molar  constant  (kmol/atm  s) 

Ko2  oxygen  valve  molar  constant  (kmol/atm  s) 

Nq  number  of  series  fuel  cells  in  the  stack 
pu2  hydrogen  partial  pressure  (atm) 

PHoO  water  partial  pressure  (atm) 

po2  oxygen  partial  pressure  (atm) 

qu?  input  molar  flow  of  hydrogen  (kmol/s) 

qo2  input  molar  flow  of  oxygen  (kmol/s) 

R  universal  gas  constant  (1  atm/kmol  K) 

Rm  stack  internal  resistance  (Q) 

T  absolute  temperature  (K) 

th2  hydrogen  time  constant  (s) 

th2o  water  time  constant  (s) 

tq2  oxygen  time  constant  (s) 


4h2 

cio. 


V 


2.2.  Reformer  model 


In  Ref.  [8],  the  author  introduced  a  simple  model  of  a  re¬ 
former  that  generates  hydrogen  through  reforming  methane. 
The  model  is  a  second-order  transfer  function.  The  mathe¬ 
matical  form  of  the  model  can  be  written  as  follows: 


<?h2  _ _ CV _ 

^methane  TiT2S2  +  (tl  +  T2)s  +  1 


where  ^methane  is  the  methane  flow  rate  (kmol/s),  CV  the 
conversion  factor  (kmol hydrogen/kmol methane),  X\,X2  the 
reformer  time  constants  (s). 

To  control  hydrogen  flow  according  to  output  power  from 
the  fuel  cell,  a  feedback  from  the  stack  current  is  considered. 
A  proportional  integral  (PI)  controller  is  used  to  control  the 
flow  rate  of  methane  in  the  reformer  [8].  Oxygen  flow  is 
determined  using  the  hydrogen-oxygen  flow  ratio  rn-o-  The 
reformer  and  the  reformer  controller  are  illustrated  in  Fig.  2, 
where  U  is  the  fuel  utilization  factor,  k 3  the  PI  gain  and  13 
the  PI  time  constant. 


2.3.  Power  conditioning  unit  model 


The  power  conditioning  unit  is  used  to  convert  dc  output 
voltage  to  ac.  As  mentioned  before,  the  power  conditioning 
unit  includes  a  dc/dc  converter  to  raise  dc  output  voltage  to 
dc  bus  voltage,  followed  by  dc/ac  inverter  to  convert  dc  bus 
voltage  to  ac.  In  this  paper,  only  a  simple  model  of  a  dc/ac 
inverter  is  considered  for  the  following  reasons:  the  dynamic 
time  constant  of  inverters  is  of  the  order  of  microseconds 
or  at  the  most  milliseconds.  The  time  constants  for  the  re¬ 
former  and  stack  are  of  the  order  of  seconds.  Hence,  includ¬ 
ing  the  inverter  time  constant  will  have  negligible  effect  on 
the  time  response  accuracy.  On  the  other  hand  it  complicates 
the  system  model.  All  commercial  FCPPs  have  to  conform 
to  IEEE  standard  no.  P-1547,  which  guarantees  the  ripples 
in  the  FCPP’s  output  voltage  to  be  within  commercially  ac¬ 
ceptable  limits.  A  simple  model  of  the  inverter  is  given  in 
Ref.  [9].  Considering  the  fuel  cell  as  a  source,  the  inverter 
and  load  connection  is  shown  in  Fig.  3.  The  output  voltage 
and  the  output  power  as  a  function  of  the  modulation  index 
and  the  phase  angle  can  be  written  as 

Vac  =  cellZ<$  (2) 


^Vcell  Vs 


sin(<$) 


Fig.  1.  PEM  fuel  cell  model. 
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Fig.  3.  Fuel  cell,  inverter  and  load  connection  diagram. 
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Q  _  (mVce ll)  -  wVceii  Vs  cos(3) 

U-  X  {  ’ 

h  = - 5 -  (5) 

Vscos(6>) 

I  =  mix  cos  (0  +  5)  (6) 

where  Vac  is  the  ac  output  voltage  of  the  inverter  (V),  m  the 
inverter  modulation  index,  8  the  phase  angle  of  the  ac  voltage 
mVceii  (rad),  Pac  the  ac  output  power  from  the  inverter  (W), 
Q  the  reactive  power  output  from  the  inverter  (VAR),  Vs 
the  load  terminal  voltage  (V),  X  the  reactance  of  the  line 
connecting  the  fuel  cell  to  the  load  (£2),  /t  the  load  current 
(A),  0  the  load  phase  angle  (rad)  and  P\  the  load  power  (W) 
(Fig.  4). 

The  dc/ac  inverter  model  is  shown  in  Fig.  4.  The  current 
feedback  signal  in  Figs.  1  and  2  is  calculated  from  the  load 
information  and  the  ac  output  voltage.  Due  to  the  assumed 
grid-independent  operation  the  transformer  model  is  not  in¬ 
cluded. 


dex  ( m ).  A  linear  activation  function  is  used  for  the  output 
layer.  The  NN  configuration  is  illustrated  in  Fig.  5. 

Two  sets  of  data  have  been  used  to  train  the  proposed 
NN.  The  two  sets  were  produced  using  the  model  given  in 
Ref.  [11].  The  following  load  profiles  were  used  to  produce 
the  training  data:  (a)  computer-simulated  step  changes  in  the 
connected  load,  (b)  actual  residential  load  profiles.  The  load 
profiles  for  a  2500  ft2  all  electric  home  with  two  adults  and 
four  children  were  used.  The  active  and  reactive  power  load 
profiles  for  a  4h  period  (4:48  p.m.  to  8:48  p.m.)  with  a  15  s 
sampling  interval  were  used. 

The  training  sets  were  used  in  a  developed  software  pro¬ 
gram  for  training  purposes  to  develop  NN  structures.  The 
developed  NN  is  in  the  form  of  weights  and  biases  to  repre¬ 
sent  the  mapping  functions  between  the  inputs  and  outputs. 
Matlab  software  is  used  to  train  the  proposed  neural  net¬ 
work  using  the  Levenberg-Marquardt  technique.  The  pro¬ 
posed  NN  and  the  FCPP  block  diagram  are  shown  in  Fig.  6. 


3.  Neural  network  controller 

In  Ref.  [11],  the  authors  introduced  a  technique  to  control 
the  active  and  reactive  power  output  from  an  FCPP.  The  pro¬ 
posed  technique  is  based  on  controlling  the  inverter  modu¬ 
lation  index  (m)  to  control  the  voltage  level  and  the  reactive 
power  output  from  the  FCPP.  The  active  power  flow  from 
the  FCPP  to  the  load  is  controlled  through  controlling  the 
phase  angle  of  the  ac  output  voltage  (5).  The  proposed  NN 
controller  uses  the  same  technique  proposed  in  Ref.  [11]  to 
control  the  active  and  reactive  power  output  from  the  FCPP. 
A  multilayer  feedforward  NN  with  back-propagation  train¬ 
ing  is  used.  The  NN  consists  of  a  fully  connected  three-layer 
network.  The  input  layer  receives  four  inputs:  from  the  load 
active  power,  load  reactive  power,  the  ac  output  voltage  mag¬ 
nitude,  and  the  hydrogen  input  to  the  FCPP.  The  hidden  layer 
has  10  neurons  with  a  log-sigmoid  activation  function.  The 
output  layer  has  two  outputs,  one  to  modify  the  ac  voltage 
phase  angle  (5),  and  the  other  to  modify  the  modulation  in- 


4.  Tests  and  results 

The  PEM  FCPP  we  currently  have  in  the  lab,  as  well  as 
the  commercial  FCPPs,  need  at  least  an  hour  to  start  from 
cold  to  build  up  the  reformer  and  stack  temperatures.  A  min¬ 
imum  load  value  (critical  load)  has  to  be  maintained  all  the 
time  during  power  plant  operation.  Decreasing  the  load  be¬ 
low  the  critical  load  value  will  put  the  FCPP  in  a  sleep  or 
dormant  mode,  where  a  small  amount  of  fuel  is  used  to  main¬ 
tain  the  reformer  and  stack  temperatures  at  operating  levels. 
Commercially  available  PEM  FCPPs  come  equipped  with 
storage  batteries  connected  in  parallel  with  the  dc  bus.  These 
batteries  serve  as  a  short-period  auxiliary  power  source  to 
meet  load  demand  that  cannot  be  met  by  the  FCPP,  par¬ 
ticularly  during  the  transient  periods.  In  the  following  test 
cases,  the  FCPP  is  assumed  to  be  in  the  active  mode  and 
the  initial  active  and  reactive  power  outputs  are  1.0  and  0.0 
respectively.  Testing  the  proposed  model  indicated  that  the 
reformer  and  the  reformer  controller  parameters  have  signif¬ 
icant  effect  on  the  power  plant  time  response.  The  reformer 
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Vac 

Pao 

Q 


parameters  affect  the  reformer  damping  factor.  Choosing  a 
damping  factor  of  1.0  for  the  uncontrolled  reformer  (open 
loop)  is  to  ensure  that  no  unrealistic  overshoot  occurs  [8]. 
This  will  yield  T\  =  T2.  For  maintaining  stability  and  rea¬ 
sonable  fuel  processor  dynamics  of  the  controlled  reformer, 
the  controller  parameters  are  chosen  to  ensure  a  damping 
factor  of  0.707  [8].  With  this  assumption  the  PI  controller 
time  constant  T3  is  equal  to  Ti  and  £3  is  equal  to  1/(2  CV). 
The  model  parameters  are  given  in  Table  1 . 

Case  1:  The  FCPP  model  shown  in  Fig.  6  is  tested  with 
step  changes  in  the  load  as  shown  in  Fig.  7.  These  abrupt 


Table  1 

Model  parameters 

Stack  temperature,  (K) 

343 

Faraday’s  constant,  F  (C/kmol) 

96484600 

Universal  gas  constant,  R  (J/kmolK) 

8314.47 

No  load  voltage,  Eo  (V) 

0.6 

Number  of  cells,  No 

88 

Kx  constant  =  No/4F  (kmol/s  A) 

0.996  x  10-6 

Utilization  factor,  U 

0.8 

Hydrogen  valve  constant,  K h0  (kmol/s  atm) 

4.22  X  10"5 

Water  valve  constant,  A^o  (kmol/s  atm) 

7.716  x  10-6 

Oxygen  valve  constant,  ATch  (kmol/s  atm) 

2.11  x  10“5 

Hydrogen  time  constant,  th0  (s) 

3.37 

Water  time  constant,  th9o  (s) 

18.418 

Oxygen  time  constant,  to2  (s) 

6.74 

Reformer  time  constant,  t\  (s) 

2 

Reformer  time  constant,  T2  (s) 

2 

Conversion  factor,  CV 

2 

Activation  voltage  constant,  B  (A-1) 

0.04777 

Activation  voltage  constant,  C  (V) 

0.0136 

Stack  internal  resistance,  Rint  (£2) 

0.00303 

Line  reactance,  X  (£2) 

0.05 

PI  gain  constants,  £5,  ko 

10 

Voltage  reference  signal,  Vr  (p.u.) 

1.0 

Methane  reference  signal,  Cmcthrcf  (kmol/s) 

0.000015 

changes  in  the  active  and  reactive  power  are  for  testing  the 
dynamic  response  of  the  system  and  do  not  necessarily  rep¬ 
resent  changes  in  a  residential  load.  The  load  model  is  cho¬ 
sen  to  reflect  all  possible  variations  of  active  and  reactive 
power.  The  change  of  current,  voltage,  output  active  and  re¬ 
active  power,  output  voltage  phase  angle,  modulation  index, 
and  flow  rate  of  hydrogen  are  illustrated  in  Figs.  8-14.  The 
time  interval  0-17  s  is  not  shown  in  the  plots  because  it  rep¬ 
resents  the  dormant  state  of  the  FCPP.  From  these  figures, 
it  is  evident  that  an  increase  in  the  load  increases  feedback 
current,  which  in  turn  decreases  output  voltage  of  the  fuel 
cell.  The  increase  in  current  increases  methane  flow  rate  and 
hydrogen  flow  rate  to  increase  phase  angle  of  output  volt¬ 
age,  which  increases  power  flow  from  the  cell  to  the  load. 
As  seen  in  the  ac  output  power  curve,  output  power  has  a 
time  delay  in  following  the  load  power.  This  is  due  to  the 
reformer  and  the  fuel  cell  time  constants.  Comparison  of  the 
results  in  Ref.  [11]  and  results  given  in  Figs.  8-14  showed 
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Fig.  7.  Load  step  changes. 
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Fig.  8.  FC  output  current. 


1.08 

1.06 

1.04 

1.02 

7  1 

a  0.98 
§  0.96 
>  0.94 
0.92 
0.9 
0.88 
0.86 


CD 

T— 

1"- 

CO 

05 

''t 

CO 

CO 

CO 

CD 

CD 

LO 

in 

"st" 

00 

CO 

o 

T— 

CM 

CO 

LO 

CD 

05 

T— 

Time  [s] 


Fig.  9.  Ac  output  voltage. 
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Fig.  10.  Active  output  power. 
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Fig.  12.  Output  voltage  phase  angle. 
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Fig.  13.  Modulation  index. 


that  the  NN  controller  was  able  to  modify  the  phase  angle 
of  the  ac  voltage  and  the  modulation  index  to  conform  to 
the  load  changes.  In  addition,  the  NN  controller  is  character¬ 
ized  by  faster  time  response  compared  to  the  PI  controllers 
used  in  Ref.  [11],  and  this  enhances  the  overall  system 
response. 

Case  2:  In  this  test  case  the  model  is  tested  using  an  ac¬ 
tual  residential  load  profile.  The  load  profile  for  the  house 
mentioned  in  Section  3  is  used.  The  active  and  reactive 
power  load  profile  for  a  2h  period  (9:00  p.m.  to  11:00 


Time  [s] 


Fig.  11.  Output  reactive  power. 


Fig.  14.  Hydrogen  flow  rate. 
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Fig.  15.  Residential  load  active  power. 


Fig.  18.  FC  reactive  power  output. 


Fig.  16.  Residential  load  reactive  power. 
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Fig.  17.  FC  active  power  output. 


p.m.)  with  a  15  s  sampling  interval  is  shown  in  Figs.  15 
and  16.  The  power  and  voltage  bases  are  5kW  and  120  V 
respectively. 

The  purpose  of  this  test  is  to  check  the  load  following 
characteristics  of  the  fuel  cell  power  plant.  The  active  and 
reactive  output  powers  are  shown  in  Figs.  17  and  18.  It  can 
be  seen  that  the  FCPP  with  the  embedded  NN  controller  has 
a  fast  response  to  residential  load  changes  and  exhibits  good 
load  following  capability. 


5.  Conclusions 

This  paper  introduces  a  technique  based  on  neural  net¬ 
work  to  control  the  active  and  reactive  power  output  from 
a  stand-alone  fuel  cell  power  plant.  The  paper  uses  an  in¬ 
tegrated  dynamic  model  for  the  fuel  cell  power  plant  and  a 
neural  network  controller  to  study  the  time  response  of  the 
fuel  cell.  The  proposed  controller  and  the  fuel  cell  power 
plant  model  have  been  tested  using  computer-simulated  step 
changes  in  the  load  active  and  reactive  power  demands,  and 
actual  active  and  reactive  load  demands  of  a  single-family 
residence.  The  results  obtained  show  rapid  response  of  the 
controller  to  load  changes,  and  the  effectiveness  of  the  pro¬ 
posed  controller  for  controlling  the  active  and  reactive  power 
output.  The  results  also  portray  the  load  following  capability 
of  the  fuel  cell  power  plant. 
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